Destruction of attractive bosonic cloud due to high spatial coherence in tight trap 
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We study coherence of trapped bosonic cloud with attractive finite-range interaction in a tight 
harmonic trap. One-body density and pair distribution function in the ground state for difi'erent trap 
sizes are calculated. We also calculate healing length and the correlation length which signify the 
presence of high spatial coherence in a very tight trap leading to the destruction of the condensate 
for a fixed particle number. This is in marked variance with the usual collapse of the attractive 
metastable condensate when N > N^r. Thus we investigate the critical frequency and critical size 
of the trap for the existence of attractive BEG. The finite range interaction gives a nonlocal effect 
in the effective many-body potential and we observe a high-density stable branch besides the known 
metastable branch. Moreover, the new branch shows universal behavior even in the very tight trap. 
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I. INTRODUCTION 

Trapped ultracold bosonic atoms have attracted a 
revival of interest since the achievement of the Bose 
Einstein condensate of dilute gases. The standard 
defining property of Bose-Einstein condensation (BEC) 
is the macroscopic occupation of the single particle 
ground state for the ideal system. In the case of 
an interacting system, BEC is often related to the 
off-diagonal long-range order [l|, By controlling the 
external magnetic field, one can virtually manipulate 
the interatomic force as well as the trap size. This can 
change the situation from a weakly interacting case to 
a very strongly interacting one and obviously the role 
of interatomic correlations come in the picture. Experi- 
mental relevance stems from the fact that experiments 
designed to produce the interference pattern in BEC 
are closely related with the study of coherence properties. 

We investigate a weakly interacting Bose gas in the 
presence of an external trap. We consider a finite 
number of bosons interacting via the finite range van 
dcr Waals interaction, instead of the commonly used 
zero-range contact interaction. The — ^ tail introduces 
an attractive non local interaction. ''Li is particularly 
interesting due to its negative scattering length and 
when the number of bosons is below some threshold 
value, a metastable condensate appears. Due to the 
attractive interaction, the effect of correlation becomes 
important even in the weakly interacting gas and we 
expect to get new physics in the study of correlation 
properties. Although it is commonly believed that the 
Gross-Pitaevskii (CP) equation, based on a mean-field 
approach and a contact interaction, is adequate for 
weakly interacting Bose gases, but in the present study 
we utilize an approximate many-body technique which 



is more rigorous and retains all possible two-body 
correlations 0-111 . This technique, called correlated 
potential harmonics (CPH) method, uses a subset of full 
hyperspherical harmonics (HH) expansion basis. While 
the latter is exact and incorporates all correlations 
in the many-body wave function, CPH basis retains 
only two-body correlations and ignores higher-body 
correlations. This assumption is manifestly valid in 
a typical laboratory BEC, which must be extremely 
dilute to avoid depletion due to recombination through 
three-body collisions. The HH method is very cumber- 
some and is practicable only for three-body systems. 
But in the CPH method, the bulk of superfluous basis 
functions are eliminated, so that it can be adopted for 
a large number of bosons (sec Sec. H for details). It 
has been successfully applied to attractive BECs and 
repulsive systems containing up to 14000 particles [1, [![. 
The method is important especially in the study of 
correlation properties in realistic condensates. Thus the 
motivation of our present work is as follows. Firstly, to 
study different correlation properties in the ground state, 
like one-body density and pair-distribution function for 
few hundred ^Li atoms in the external three-dimensional 
trap. This facilitates the investigation of the effect of 
interatomic correlations in the many-body wave function 
on the bulk correlation properties of the condensate. To 
quantify the latter, we also calculate the healing length 
which is another key quantity in the study of coherence 
properties. Secondly, we investigate the effect of the 
change of trap size on the correlation properties. By 
controlling the trap frequency it is possible to make 
tighter traps. We observe that in a tighter trap the 
metastable condensate becomes highly correlated. Ef- 
fective correlation length and the healing length reduce 
drastically, as the trap size decreases; however, they 
remain larger than the trap size, indicating strong bulk 
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spatial correlation within the metastable condensate. In 
addition to the metastable condensate, the many-body 
system exhibits a strongly attractive narrow well, which 
supports cluster states (see below). Below a critical trap 
size the metastable region in the effective many-body 
potential disappears. This manifests the destruction of a 
metastable condensate in a very tight trap. Wc estimate 
the critical trap size above which BEC can be observed. 
The disappearance of the condensate in a very tight trap 
significantly differs from the usual collapse of attractive 
BEC in a large trap. 



Due to the presence of a realistic interaction with 
a strong repulsive core, we have a deep but finite 
attractive well outside a repulsive core on the left side of 
the metastable region in the effective potential. Hence, 
besides the appearance of the low density metastable 
BEC, we also find a stable branch at high density 
which leads to the formation of atomic clusters. With 
decreasing trap size, the metastable branch shrinks and 
eventually disappears (resulting in a collapse due to 
squeezing). However the attractive well of the effective 
potential supporting the high density branch remains 
invariant in position, magnitude and shape with decrease 
in trap size. Thus it exhibits a universal behavior. The 
energy and size of the cluster do not change with the 
trap size. In earlier calculations j^-Q, a non-local 
effective interaction in the momentum space, which is 
equivalent to a local repulsive contact interaction plus 
a finite-range attractive interaction in the coordinate 
space was adopted in a variational calculation of the 
GP equation. This also produced the narrow attractive 
well and the associated high-density branch, in addition 
to the low-density metastable branch. On the other 
hand, a purely attractive local contact interaction 
in the GP equation produces a pathological essential 
singularity at the origin, besides the metastable region. 
In Ref . [y, 0] , quantum tunneling rates between the low- 
and high-density branches and the decay rates due to 
two- and three-body collisions were reported. Nature 
of transitions between the two branches as a function 
of the number of bosons was also investigated. In the 
present study, we obtain these two branches as a result 
of using the realistic van dcr Waals interaction having a 
strong short-range repulsion (due to nucleus-nucleus re- 
pulsion) and a long-range attraction, in an approximate 
many-body theory, which incorporates all two-body 
correlations. We focus on the correlation properties and 
investigate the effect of decreasing trap size. 



The paper is organized as follows. The many-body 
calculation technique is presented in Sec. II. Sec. HI 
presents the calculation of many-body effective poten- 
tial and different correlation properties like one-body 
density and pair distribution function. Sec. IV contains 
a summary and the conclusions. 



II. MANY-BODY CALCULATION WITH 
CORRELATED POTENTIAL HARMONICS 
BASIS 

The Hamiltonian for a system of {N + 1) atoms ( each 
of mass m) and interacting via two-body potential has 
the form 



■ 2 N+1 N+1 
i—\ z>j — 1 



JV-l-1 



i>j=l 



. (1) 

where V{xi — Xj) = Vi^ij) is the two-body potential de- 
scribed later and Xi is the position vector of the i-th 
particle. We use the standard Jacobi coordinates defined 



as Ci 



[i = l,2,....iV) and 
TiT.f=iX,. Then the 



1 -* 
Xi+1 ^ - l^j=l Xj 

the center of mass through R — 

relative motion of the atoms is described in terms of N 
Jacobi vectors (Ci, ■ ■ Ov) ^ [1| 

N 

■ ^ V^, + Vtrap{r) + V{Ci, ■ ■ •, Ov) - E 



fi2 

m 



xvI-(Ci,...,Cw)==0, (2) 

where V is the sum of all pair- wise interactions expressed 
in terms of the Jacobi vectors. We assume that only two- 
body correlations in the many-body wave function are 
important 0, 0] , This permits the use of the potential 
harmonics expansion method, in which the total wave- 
function 4* is decomposed into two-body Faddeev com- 
ponents tpij for the (ij) interacting pair. 



N + 1 
ij>i 



(3) 



Note that with two-body correlations alone, ^pij is a 
function only of two-body separation (r^ ) and a global 

length, called hypcrradius {r = Cf]) a-nd is inde- 

pendent of the coordinates of all the particles other than 
the interacting pair ipij (symmetric under Pij) 

satisfies the Faddeev equation 

-V{n,) Y ^ki{x), (4) 

k,l>k 



(T -I- Vtrap - E)4ly{x) 



T being the total kinetic energy; operating ^-^ • on 
both sides of Eq. wc get back the original Schrodingcr 
equation. In this approach we assume that when [ij) 
pair interacts, the rest of the bosons are inert spec- 
tators. Moreover, since particle labels are unimpor- 
tant, we take as Qn- Next we define a hypcrradius 



Pij 



'Ylik=i Cfc for the remaining (A^ — 1) noninter- 



acting bosons such that 



Pi 



and Tij = r cos ( 
dinates are 



Ptj 



r smt 



In this choice the hyperspherical coor- 



(r, 17 at) = (r, (/), 1?, ip, VtN-i) 



(5) 
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where ip) are the two spherical polar angles of the 
separation vector r^, n^v-i involves (3A^ — 4) variables: 
2{N — 1) polar angles associated with [N — 1) Jacobi vec- 
tors Cii ■ • Cjv-1 and [N — 2) angles defining the relative 
lengths of these Jacobi vectors [31 ■ Then the Laplacian 
in 3A^-diniensional space has the form 



N 



2,N -1 d ^ L^inpf) 
r dr J.2 ■ 



(6) 



L^(rijv) is the grand orbital operator in D = 3A^ dimen- 
sional space. Eigcnfunctions of this operator, correspond- 
ing to all possible sets of quantum numbers (associated 
with 3A^— 1 degrees of freedom), constitute the complete 
hyperspherical harmonics (HH) basis. Expansion of ^ 
in this basis would lead to an exact treatment. How- 
ever, the degeneracy of the basis (arising from different 
allowed values of iN — 1 quantum numbers for a particu- 
lar grand orbital quantum number, C) increases very fast 
with C and N . This fact and the fact that calculation of 
matrix elements involve (3A'^ — l)-dimensional integrals 
make the calculation extremely cumbersome for iV > 3. 
Moreover, imposition of symmetry also becomes very dif- 
ficult as N increases. Hence the HH expansion method 
(HHEM) is practicable only for three-body systems. A 
great deal of simplification is possible for the laboratory 
BEC, which is extremely dilute and only two-body corre- 
lations are important. One can then use a subset [called 
potential harmonics (PH) basis] of HH basis with great 
advantage. Potential harmonics for the (ij)-partition are 
defined as the eigenfunctions of L'^{Q,m) corresponding to 
zero eigenvalue of L^I^Im-i) ■ The corresponding eigen- 
value equation satisfied by L'^{VIn) is Q 

[L'i^N) +C{C + D- 2)] r^iK^i{n,j) = 0, C = 2K+1 

(7) 

This new basis is a subset of the full HH set and it does 
not contain any function of the coordinateCi with i < N. 
It is given by a simple closed expression [9| 

'P2K+li^iv)) - YlmM ^''^P^K+MyoiD - 3), (8) 

where Yim{ujij) is the spherical harmonics and Uij ~ 
("!?, (/s). The function '^^^ P^^k^tM') expressed in terms 
of Jacobi polynomials [l^ and J^o (3 — 3) is the HH of 
order zero in the (3iV — 3) dimensional space, spanned by 

{Ci) • • • ; Cjv-i} Jacobi vectorsjl]. Thus the contribution 
to the grand orbital quantum number comes only from 
the interacting pair and the 3A dimensional Schrodingcr 
equation reduces effectively to a four dimensional equa- 
tion. The relevant set of quantum numbers (associated 
with the hyperangles) are only three - orbital Z, azimuthal 
771 and grand orbital 2K -I- I for any N . This leads to a 
dramatic simplification of the many-body calculations. 
Besides drastic reduction of degeneracy of the basis, po- 
tential matrix elements involve only three-dimensional 
(one-dimensional for central forces) integrals. The phys- 
ical picture is that all irrelevant degrees of freedom have 



been frozen. Using this procedure, we have solved BEC 
containing up to 14000 bosons d, . The method has 
been successfully applied to attractive condensates as 
weh [5l. [l^ - [l^ . We expand (ij)-Faddeev component, 
in the complete set of potential harmonics appropriate for 
the (ij) partition: 



E-plva (Q{if. 



''K 



(r). 



(9) 



K 



Note that the notation has been slightly modified to in- 
clude the superscript [ij) in ilTv, to indicate that it is 
the full set of hyperangles in D dimensional space, for 
the particular choice of Jacobi vector — ^ij. Eq. ([9]) 
includes two-body correlations only. This is perfectly jus- 
tified in the context of dilute attractive Bose gas, where 
the effect of two-body correlation is important and one 
can safely ignore the effects of higher-body correlations. 
Taking projection of Eq. on a particular PH, a set of 
coupled differential equations (CDE) is obtained [3, Q 

TT + + 1) + ^K{K + a + p + l)] 

777, dr^ 77ir^ 

+ Vtrap{r)-E\uKi{r) (10) 



Y,fKlVKK'{r)fK'lUK'l{r) 



K' 



3Ar-3 



3Ar-5 



where UKi{r) = fKiu'-^^ir), C = 1+ _ 
/3 = I + ^, I being the orbital angular momentum con- 
tributed by the interacting pair and K is the hyperangu- 
lar momentum quantum number. The potential matrix 
element Vkk' {r) is given by [1, El 

VKK'ir) = [ ^^'^;,(f^^^'^)^(x.,)7'^'?J,+,(l^^^Vf^^''^- 



(11) 

The quantity is given by 

fh = E < ^2/?+;(f^!;'')i^2/?+;(<^) > • (12) 

k,l>k 

It is the overlap of the PH for the (ij)-partition (corre- 
sponding to only the (ij)-pair interacting) with the sum 
of PHs for all partitions. An expression for /|-; can be 
found in Ref Q. So far we have disregarded the effect 
of the short range correlation in the PH basis. In the 
case of a dilute Bose gas, as the energy of the interact- 
ing pair is extremely small, the two-body interaction is 
generally represented by the s-wave scattering length a^c 
alone (as in the mean- field GP treatment), disregarding 
its detailed structure. On the other hand, a realistic in- 
teratomic potential like the van der Waals interaction has 
an attractive — \- tail at larger separations and a strong 

repulsion at short separation. The short-range behavior 
is usually represented by a hard core of radius r^- For a 
given two-body potential having a finite range, a^c can be 
obtained from the asymptotic solution of the zero-energy 
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two-body Schrodingcr equation 



h^J^ d_ 

m rfj drij 



dri{ri 



+ VinMrv) = 0. (13) 



The value of rc is obtained from the requirement that 
the calculated asc has the experimental value. The 
correlation function quickly attains its asymptotic form 
(Ci -|- C2/rij) for large r^. The asymptotic normaliza- 
tion is chosen to make the wavefunction positive at large 
rij and the corresponding scattering length is given by 
asc = ~§7- In the experimental BEC, the energy of the 
interacting pair is negligible compared with the depth 
of the interaction potential. Thus 'r]{rij) correctly repro- 
duces the short rij behavior of ipij (rij , r) . Hence we in- 
troduce this as a short-range correlation function in the 
expansion basis and call it as correlated potential har- 
monics (CPH) basis Q. 



2K+1 



(f](,,))j =y,„(^,,)^'^^^2K+/W 

-I cori^ 

xyo{3N-3)7j{n,), (14) 

The correlated potential matrix element Vkk' (r) is now 
given by 



1 + z 



xP^f(z),y 




(15) 



Here h'^ and wi{z) are respectively the norm and weight 
function of the Jacobi polynomial P^^{z) [1,01 . 
Note that the inclusion of the short-range correlation 
function, r](rij) makes the PH basis non-orthogonal. 
This introduces an overlap matrix on the eigenvalue side 
of the matrix eigen value equation. One can use the 
standard procedure for handling a non-orthogonal ba- 
sis, by introducing a transformation using the eigen- 
values of the overlap matrix to convert the eigenvalue 
equation into the standard diagonalization of a sym- 
metric matrix. However, we have checked that r]{rij) 
differs from a constant value only by a small amount 
in a relatively small interval. Hence to simplify the 
calculation, we project Eq. (4) on to a particular PH, 

viz. 'P2K+i^^%^'^)- The dependence of the overlap < 

Vi'^+i{^N^)\Vi'^c+ii^f^)v[rki) > on the hyperradius r 
is quite small. Disregarding derivatives of this overlap 
with respect to the hyperradius, we approximately get 
back Eq. (10), with VKK'{r) given by Eq. (15). 



III. RESULTS 

A: CHOICE OF INTERACTION AND 
CALCULATION OF EFFECTIVE POTENTIAL 

We consider the ^Li condensate with asc = —27.3 
Bohr. For the numerical calculation, we choose the oscil- 



lator units (o.u.) of length and energy, commonly used 
in BEC calculations. The oscillator length, aho ~ 
is chosen as the unit of length and the oscillator energy 
hcu as the unit of energy, uj being the trapping frequency 
(w = 2t:i'). However for presenting results, both in tab- 
ular form and in figures, we use MKS units: meter (m) 
and Joule (J) for length and energy respectively. The 
trap size corresponding to Rice University experiment is 
o-ho = 3.0/im. Our chosen potential is the realistic van der 
Waals potential, which has a strong repulsive core (which 
is chosen as a hard core of radius Tc) and an attractive tail 
at larger separations: V{rij) = oo for rij < rc and = — 

for rij > rc- The strength (Cg) is known for a given type 
of atom. In the limit of Ce — >■ 0, the potential becomes a 
hard sphere and rc coincides with the s-wave scattering 
length Use- For the potential including the long range 
part, a tiny change in rv may cause an enormous change 
in Use, including sign [5|. As Tc decreases from a large 
value, Use decreases, and at a particular critical value of 
Tc, it passes through an infinite discontinuity from — oo 
to -f oo 0. Thereafter the potential supports a two body 
bound state. This pattern repeats as rc decreases further. 
Positive values of asc correspond to repulsive potential 
whereas negative Usc values correspond to attractive po- 
tential. Thus minute tuning of rc can cause the effective 
potential to change from attractive to repulsive. In the 
mean-field description, the two-body interaction is solely 
represented by the s-wave scattering length asc and de- 
pending on its sign, positive or negative, the condensate 
is treated as repulsive or attractive respectively. How- 
ever in our many-body calculation, we solve zero-energy 
two-body Schrodinger equation with V{rij) given above 
and tune rc to obtain correct value of a^c which mimics 
the ''Li condensate of Rice University [l^ . However with 
a tiny change in r^ asc may change by large amount in- 
cluding the sign. Each additional sign change means that 
the potential will support an extra two-body bound state 
and it results in an extra node in ij{rij). We choose rc, 
such that it corresponds to the zero node in the two-body 
wave function. The chosen parameter for our calculation 
is Ce = 1.71487 x IQ-^^ o.u. and rc = 5.3378 x 10"'' o.u.. 
With these sets of parameters we solve the set of coupled 
differential equations by the hyperspherical adiabatic ap- 
proximation (HAA) In HAA, one assumes that the 
hyperradial motion is slow compared to the hyperangu- 
lar motion. The effective potential for the hyperradial 
motion (obtained by diagonalizing the potential matrix 
together with the diagonal hypercentrifugal repulsion for 
each value of r) is obtained as a parametric function of r. 
We choose the lowest eigenpotential (a;o(r)) as the effec- 
tive potential. Thus in HAA, energy and wavefunction 
are obtained approximately by solving a single uncoupled 
differential equation 



?l2 ^2 

+^oW 

m dr^ 



E 

K=0 



dr 



E 



Co(r) =0, 
(16) 
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subject to appropriate boundary conditions on Co('")- 
The third term is a correction to the lowest order HAA 
approximation. XKoi^) is the A'-th component {K 
being the hyperangular momentum quantum number) of 
the eigenvector corresponding to the lowest eigenvalue 
{LOo{r)) of potential plus hypercentrifugal matrix. This 
is called uncoupled adiabatic approximation (UAA), 
whereas disregarding the third term corresponds to the 
extreme adiabatic approximation (EAA). The principal 
advantages of the present method are as follows. First, 
the correlated PH basis keeps all possible two-body cor- 
relations and the number of variables is reduced to only 
four irrespective of the number of bosons in the trap. 
So by this method, we can treat quite a large number of 
atoms without much computational difficulty. Second, 
the use of HAA basically reduces the multidimensional 
problem to an effective one-dimensional one introducing 
the effective potential. The effective potential (wq) gives 
clear qualitative as well as quantitative pictures. For 
our numerical calculation we fix / = and truncate the 
CPH basis to a maximum value K = Kmax, requiring 
proper convergence. Third, the use of van der Waals 
interaction with a finite range is more realistic than 
the use of a zero-range contact interaction, as in the 
mean-field GP approach. The pathological singularity 
of the (5-function attractive potential does not arise in 
the present treatment. 

In Fig. 1, we plot the many-body effective potential 
ojQ^r) as a function of hyperradius r for A^=200 atoms in 
the condensate. As iV = 200 is less than Ncr ^ 1300 fis} . 
the condensate is metastable and is associated with a 
deep and narrow attractive well (NAW) on the left side. 
For r — > 0, there is a strong repulsive wall. This is the 
immediate refiection of using hard core van der Waals 
interaction. The effective potential for an attractive 
contact interaction goes rapidly to — oo as r — > 0, 
causing an essential and pathological singularity. Thus, 
our many-body picture with nonlocal interaction is in 
sharp contrast with the GP mean-field picture with 
local interaction. The nonlocal interaction and the 
repulsive core of the van der Waals potential prevent 
the Hamiltonian from being unbound from below. For 
N less than the critical value Ncr, a metastable region 
(MSR) appears for larger r, an intermediate barrier (IB) 
separating the NAW and the MSR. For still larger r, 
the influence of the attractive interaction subsides and 
the external wall of the harmonic trap dominates. In 
panel (a) of Fig. 1, the NAW together with the repulsive 
core is shown. In panel (b) of the same figure, the IB 
and MSR have been plotted. The bottom of the NAW 
has a very large magnitude compared with the bottom 
of the MSR, hence they cannot be shown in the same 
figure. Note the widely different scales used in the two 
panels. Furthermore, r in panel (b) is in logarithmic 
scale. 

With the increase in the number N, we observe 



a decrease in the height of the intermediate barrier, 
together with a decrease in the difference between the 
maximum of IB and the minimum of MSR and the NAW 
starts to be more negative and narrower. As iV — >■ Ncr, 
the maximum of IB and the minimum of MSR merge to 
form a point of inflexion, with the disappearance of the 
MSR. At = Ncr, the metastable condensate collapses. 
For N > Ncr, there will be only the NAW and no 
metastable condensate. In our present study we observe 
the IB just vanishes and the condensate collapses at 
N = 1430. This is the usual collapse of the attractive 
condensate. At this point, all the atoms get trapped 
in the NAW and form van der Waals cluster which 
corresponds to a high-density branch in the density 
profile. 

Next we study different correlation properties in such 
a realistic condensate. We also investigate how the 
coherence properties depend on the trap size (a/jo) and 
determine the critical size of the trap. 



B : ONE-BODY DENSITY 

It was originally pointed out that the BEG is evidenced 
by the presence of off-diagonal long-range order in the 
one-particle density matrix [l|. However this definition 
is not strictly valid in the case of a finite-size inhomoge- 
neous system. Due to the presence of an external har- 
monic trap, our system is inhomogeneous. In our many- 
body formalism, we define the one-body density as the 
probability density of finding a boson at a distance rk 
from the center of mass of the condensate as 

Ri{rk)^ jji^^dr'. (17) 

where ^ is the full many-body wave function and the 
integral over the hypervolume r' excludes the variable 
Tfc. The incremental hypervolume dr' is given by 

dr' = r'^^-^ cos^ ^sin^^^^ (jjdr' d(l)duj^jd^N-2 (18) 

where r' is obtained from the relation 

^ r'^ + 2rl (19) 

and the other symbols have their usual meanings [iTj . 
The integral is computed analytically followed by numer- 
ical computation using a 32-point Gaussian quadrature 
with the original interval divided into progressively 
increasing subintervals. According to Penrose and 
Onsager definition of Bose-Einstein condensation, the 
one-particle density matrix must be associated with a 
single macroscopic eigenvalue [l|. However the function 
ipif) can be directly expressed in terms of first-order 
correlation function [l^ and all the information of the 
one-body density correlation are contained in Eq. (17). 
The short-range repulsion between interacting atoms 
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give some new aspects in the atomic correlation. We 
present our results in Fig. 2 for 200 ^Li atoms with 
flsc = —27.3 Bohr and two different trap sizes: trap 
frequency uj = 1.01 kHz (corresponding to trap size 
o-ho = 3.0 /im) and lj = 2.27 kHz (corresponding to 
trap size ato = 2.0 /.jm). The one-body density profile 
deviates from the Gaussian profile of non-interacting 
case. For comparison, we include the mean-field GP 
results. The deviation from the GP result is attributed 
to the effect of interatomic correlation. For the re- 
duced trap size of 2.0 ^m, we observe that the peak 
becomes higher and narrower, the difference between 
the mean-field GP and many-body result also increases, 
as the system develops shorter range correlations. This 
can be understood from the fact that in the standard 
GP approach the NAW is disregarded, which pulls the 
system inwards. The many-body results for gradually 
reduced trap sizes are presented in Fig. 3. It was 
observed earlier that for the homogeneous system the 
long-range behavior in the one-body density follows a 
power-law decay [Tsj . However in the presence of an 
external trap the long-range tail in the one-body density 
is of the order of the trap width. Thus by reducing 
the trap size gradually we observe that the long tail 
in the one-body density decreases as expected. For a 
very tight trap, the one-body density is sharply peaked 
and the long-range order is sharply reduced, in tune 
with reduced trap size. It indicates that the atoms in 
the metastable condensate are highly correlated. These 
features become more clear from Fig. 4, where we plot 
the change in the MSR of the effective many-body 
potential with the change in the trap size. In Table I 
we present the position and the value of the maximum 
of IB and the second minimum (which comprises the 
MSR) in the column 3 and 4 respectively. Although 
the position and the value of the maximum of IB do 
not change much with the trap size, the position and 
the value of second minimum are greatly shifted. It 
indicates that the MSR is pulled in and the condensate 
shrinks with decrease in trap size. The corresponding 
interaction energy is presented in column 5 of Table I. 
With decrease in trap size the attractive interaction 
energy sharply increases and we find auo = 0.42 ^m 
as the critical trap size. Just below this critical size 
the condensate will be destroyed due to high quantum 
fluctuations as the negative interaction energy increases 
very fast for smaller auo- What actually happens is the 
following. As the trap size is reduced, the system gets 
squeezed. Consequently, the interatomic spacing reduces 
and the net attractive interaction energy increases. 
Eventually, when the interatomic spacing reaches the 
typical cluster size, tightly bound clusters are formed, 
together with the complete removal of the metastable 
condensate. Such a scenario is not possible with an 
attractive contact interaction, for which the effective 
potential goes to — oo as r — >■ 0. Note that this collapse 
is different from the usual collapse for a^o = 3/j,m, where 
the metastable region vanishes at iV = Ncr ■ The critical 



trap size strongly depends on the number of bosons in 
the trap, which may be just few tens to a few hundreds 
depending on the trap frequency. For N = 200 atoms 
the critical trap frequency is u = 51.47 kHz. 

The healing length is sometimes referred to as the 
coherence length and may be considered a relevant 
quantity to quantify the correlation in very tight traps. 
Healing length ^ is basically the minimum distance over 
which the order parameter can heal. It is in general 
calculated by balancing the quantum pressure and the 
interaction energy of the condensate. In Fig. 5 we plot 
the healing length as a function of the trap size, which 
shows a steep decrease in Q in very tight traps. The 
simple expression for healing length is ^ = . , 

^ 8irascn(0) 

where n(0) is the peak density at r = 0. However 
in principle n depends on r and Cir) ~ . ^ , 

y/ 87Tascn{r) 

i.e., as n{r) 0. (^(r) — > oo at the surface of the 
cloud. We can also compare our results with the simple 

J, 2 

Gaussian estimate, for which n{r) = — . ^ e "^o and 

C(0) ~ al^^^ s/Sira^cN . Thus C scales as a,' . Thus the 
smooth decrease of the healing length in Fig. 5 reflects 
the effect of interatomic interaction and two-body 
correlation. 

As mentioned earlier, in the standard GP mean-field 
theory, the metastable region vanishes as iV — ?> Nc- and 
the condensate collapses into the singular well. The fate 
of the condensate is not predicted further. Here, the 
many-body picture is different. In our case the narrow 
attractive well on the left facilitates further study of 
the condensate near the criticality. The metastable 
condensate leaks through the intermediate barrier, and 
settles down in the NAW. The atoms in the condensate 
form van der Walls clusters via strongly enhanced 
three- and higher-body collisions. It corresponds to 
the high-density branch jToj . In Fig. 6 we demonstrate 
graphically the dependence of the NAW on the trap 
size. In Table I, we present the position and value 
of the deep minimum of the NAW in column 1. The 
position of the deep well and the adjacent barrier, 
which comprise the NAW, do not change with the trap 
size. These clearly demonstrate that the high-density 
branch remains invariant with the trap size. The calcu- 
lated size of the atomic cluster is of the order of 0.005 /im. 

This high-density stable state within the NAW cor- 
responds to the Li clusters. It would be interesting 
if these clusters could be detected experimentally. Al- 
though there is no experimental study of "^Li^ clusters 
till now, bosonic '^Hcn clusters have been studied ex- 
perimentally. The size of ^He clusters is determined 
by the matter wave diffraction gratings [2^. This is 
a very promising technique for the study of ^HeAr clus- 
ters. However, such techniques may not be feasible for 
clusters formed from collapsed BEG, due to small num- 
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ber of such clusters and difBculty in extracting a beam 
for the diffraction experiment. Perhaps a better method 
may be using spectroscopic techniques [llj . Utihzing the 
Feshbach resonance the effective interaction between two 
atoms can be changed essentiaUy to any value as desired 
and it facilitates the creation of large weakly bound clus- 
ters. The signature of the universal behavior of weakly 
bound bosonic cluster is studied in Ref. [22| and can be 
observed in ultracold Bose gas. 



tighter. During this squeezing process, the correlation 
length remains larger than the trap size. It indicates that 
the metastable condensate remains highly correlated in 
the tight confinement. But when it is squeezed beyond 
the critical trap size, clustering of particles occurs with 
the destruction of metastable BEC. 



IV. CONCLUSION 



C: PAIR-DISTRIBUTION FUNCTION 

So far we have focused on the one-particle density 
which basically considers analysis of the order param- 
eter and contains all information about the one-particle 
aspects. In the present section we are interested to calcu- 
late the probability R2 {rij ) which is defined as the prob- 
ability of finding the (ij)-pair of particles at a relative 
separation rij. For a correlated interacting system we 
define R2{rij) as follows 

R^inj) = \yj\'dr", (20) 

where ip is the many-body wavefunction as before but 
now the integral over the hypcrvolumc r" excludes in- 
tegration over Tij. The incremental hypcrvolumc dr" is 
given by 

dr" = rfj{rs\n(j)f^-^duJijdpijdnN-i (21) 

Our results for the same number of ^Li atoms and for the 
same trap sizes as in Fig. 3 are presented in Fig. 7. When 
Rai^ij) =0, there is no diagonal correlation. The peak 
value of i?2 at some distance signifies strong cluster- 
ing effects. The study of pair correlation is important 
as the realistic interatomic interaction with a small hard 
core repulsion plays a crucial role. It forbids the atoms 
to come too close due to nucleus-nucleus repulsion. Thus 
R2{fij) is always zero at = 0. Unlike the uniform 
system with no confinement, our calculation shows that 
R2 vanishes asymptotically. This is the effect of exter- 
nal confinement which restricts the pair separation to 
a finite value. In our earlier calculations we observed 
the dependence of correlation function on the interaction 
strength However here we observe that even for 

weak interaction and just a few hundred atoms, correla- 
tion length decreases drastically in a tight trap. To be 
more quantitative, we define the correlation length as the 
width of correlation function at the half of its maximum 
value and plot it as a function of trap size in Fig. 8. In 
Fig. 8 we also plot the position of clustering spot as a 
function of trap size. Clustering spot is defined as the 
position of where pair-distribution function attains 
a maximum value. Even for such weak interaction, we 
observe that both the correlation length and position of 
cluster spot in the matastable condensate decrease dras- 
tically, in tune with the trap size, as the trap becomes 



In an attractive BEC, the destruction of the conden- 
sate takes place when the number of bosons in the exter- 
nal trap exceeds the critical number. At this point, the 
interaction energy becomes large negative (so that the 
kinetic pressure fails to balance it) and the condensate 
shrinks, such that the density at the center of the trap 
becomes very high. The ultimate fate is the collapse of 
the attractive BEC. As the system becomes highly corre- 
lated, it is appropriate to undertake a correlated many- 
body approach to describe several correlation properties 
of a realistic condensate in a harmonic confinement. The 
use of a realistic interatomic interaction and the finite 
size trap give the realistic features in one-body density 
and pair correlation properties. We calculate such corre- 
lation properties of an attractive BEC, using the corre- 
lated potential harmonics technique. As the trap size is 

given by a^o ~ \J'^^i it is easy to change the trap size 
by controlling the laser frequency. A high trapping fre- 
quency leads to the creation of a very tight trap. In this 
report, we undertake the study of correlation properties 
in a tight trap. Our realistic approach produces an effec- 
tive potential in which the condensate moves. This effec- 
tive potential for a metastable condensate has a strongly 
repulsive core, followed successively by a narrow attrac- 
tive well, an intermediate barrier, a metastable region in 
which the metastable condensate resides and finally the 
high outer wall of the trap. This is in sharp contrast with 
the mean-field CP approach, for which there is an attrac- 
tive singularity on the left of the intermediate barrier. 
The present approach produces a realistic post-collapse 
scenario. We observe that in a very tight trap, the sys- 
tem becomes highly correlated and effective correlation 
length is drastically reduced. Even for a small number 
of atoms well below the critical number, the metastable 
condensate collapses when the trap size is below a certain 
value. One can define a critical trap size and critical trap 
frequency for the existence of a metastable BEC in tight 
confinement for a given number of atoms with attractive 
interaction. Below the critical size, the metastable con- 
densate will be destroyed due to high spatial coherence 
and all the atoms will settle down in the narrow attrac- 
tive well, with the formation of coherent van der Waals 
clusters. For small trap sizes above the critical value, 
the narrow attractive well shows a universal character, 
becoming independent of the trap size. This shows that 
the clusters cannot be further squeezed after their forma- 
tion. 
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FIG. 1. (Color online) Plot of the effective potential ioo{r) 
against r for 200 ^Li atoms in the usual trap of size a^o = 
3.0 fim. The upper panel shows the narrow attractive well 
(NAW). The lower panel shows the metastable region sepa- 
rated from NAW by the intermediate barrier (IB). Note the 
different scales used for r and CL)o(r) in the two panels, so that 
NAW is far to the left of the plot in panel (b) in which r 
is in logarithmic scale. Therefore, the NAW is not visible in 
panel (b). Note also that the bottom of the NAW would have 
a value ~ —3 x 10* in the unit used in panel (b). 



TABLE I. Parameters for NAW, IB and MSR of the effective 
potential and the interaction energy < V > for different trap 
sizes. Ti and uji correspond to the position and value of the 
extrema of NAW (i = 1), IB [i = 2) and MSR (i = 3). Note 
that different units have been used in different columns. 
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NAW 


IB 


MSR 


<V > 












(Aim) 


(10-**m) (IQ-^^J) 


(Aim) (10-=*J) 


{fj.m) (10-=»J) 


{10-'^" J) 


3.0 


5.79 -2.93 


5.40 9.56 


71.70 0.31 


-0.861 


2.0 


5.78 -2.92 


5.40 9.54 


47.40 0.68 


-2.915 


1.0 


5.80 -2.92 


5.50 9.61 


22.80 2.61 


-26. 127 


0.5 


5.80 -2.92 


5.70 10.72 


10.05 9.22 


-306.84 


0.42 


5.80 -2.92 


6.47 12.11 


7.22 12.09 


-1321.94 



EC 
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0.015 
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a^o = 3.0 (im (GP) 
aho = 2.0 m (PH) 
aho = 2.0 m (GP) 



200 Li atoms 
age = - 27.3 Bohr 




FIG. 2. (Color online) Plot of one-body density (i?i(rfc)) 
against for 200 ^Li atoms in the trap with the size of 3.0 
fim and 2.0 fim. Comparison with the mean field GP results 
is also presented. 
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FIG. 3. (Color online) Plot of one-body density (Ri{rk)) 
against Vk for the same number of atoms as before with dif- 
ferent trap sizes. 
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FIG. 4. (Color online) Change in the position and the depth 
of metastable region for 200 atoms with different trap sizes. 
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FIG. 5. (Color online) Plot of healing length as a function of 
trap size. 
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FIG. 6. (Color online) Plot of narrow attractive well for var- 
ious trap sizes. 




FIG. 7. (Color online) Plot of pair distribution function 
{R2{rij)) against Vij for an attractive Bose gas with differ- 
ent trap sizes. 




FIG. 8. (Color online) Plot of correlation length and cluster- 
ing spot as a function of trap size. 



